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ABSTRACT 


The  set  of  hydrodynamic  equations  together  with  Maxwell's 

equations  which  describe  a  nonuniform  plasma  as  formulated  by  Cohen  J19B2J 

have  been  applied  to  a  two  dimensional  problem  where  the  magnetic  field 

has  only  one  component  which  is  transverse  to  the  direction  of  propaga- 

-Bx 

tion  and  the  profile  is  of  the  form  e  .  These  assumptions  lead  to  a 
set  of  second  order  coupled  electromagnetic  and  hydrodvnamic  differen¬ 
tial  equations  which  describe  the  region  under  consideration.  Solutions 
to  this  coupled  set  have  been  obtained  by  the  application  of  the  method 
of  Frobenius.  Considerable  simplif ication  resulted  in  the  solutions 
when  the  assumption  that  the  rms  velocity  of  the  electrons  is  much 
smaller  than  the  speed  of  light  was  made.  The  resulting  set  of  approxi¬ 
mate  solutions  have  been  applied  to  the  problem  of  obtaining  reflection 
and  transmission  coefficients  for  the  case  of  a  plane  wave  incident  from 
free  space  upon  a  layer  of  compressible  inhomogeneous  plasma. 
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1.  INTRODUCTION 


While  problems  dealing  with  compressible  homogeneous  plasmas 
and  incompressible  homogeneous  and  inhomogeneous  plasmas  have  receiveo 
considerable  attention  in  the  past,  little  has  appeared  on  the  problem 
of  propagation  through  a  plasma  that  is  both  compressible  and  inhomo¬ 
geneous.  Furthermore,  one  finus  that  frequent  use  is  made  of  the  bound¬ 
ary  condition  that  tne  normal  component  of  tne  electron  velocity  shall 
vanish  at  the  free  space-plasma  interface.  While  the  application  of 
this  boundary  condition  may  yield  useable  answers,  it  is  doubtful  tnat 
rigid  interfaces  can  be  realized  in  situations  dealing  with  reentry 
communications.  Tne  boundary  condition  in  question  can  be  avoided  by 
suitably  choosing  a  profile  which  is  continuous  at  the  air  plasma  inter¬ 
face.  The  electron  density  will  consequently  be  zero  there  and  tne 
vanishing  of  the  normal  component  of  tne  electron  velocity  will  be  sat¬ 
isfied  automatically.  Additional  boundary  conditions  for  the  pressure 
and  its  first  derivative,  however,  must  be  introduced  in  order  to  specify 
tue  problem  uuiquely. 

In  the  following,  the  equations  describing  a  compressible  in¬ 
homogeneous  plasma  will  be  derived  first.  Then  a  special  case  of  tne 
latter  will  be  considered  and  the  corresponding  solutions  will  be  found. 
Lastly,  these  will  be  applied  to  the  determination  of  the  reflection  and 
transmission  coefficients  for  the  case  of  a  plane  wave  incident  upon  a 
layer  of  compressible  inhomogeneous  plasma. 
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2.  STATEMENT  OF  THE  PROBLEM 

Let  us  focus  our  attention  on  the  problem  of  determining  the 
transmission  and  reflection  coefficients  for  the  case  of  a  plane  wave 
incident  obliquely  upon  a  plasma  layer  which  is  both  inhomogeneous  and 
compressible.  Let  the  layer  be  bounded  by  free  space  for  the  region 
x<0  on  one  side,  and  on  the  other,  for  x>t,  by  a  compressible  homo¬ 
geneous  plasma.  The  profile  describing  the  inhoaogeneity  of  the  layer 
will  be  of  the  form  exp(-Bx).  The  relative  dielectric  constant  of  free 
space  will  be  taken  as  unity  and  that  of  the  compressible  homogeneous 
region  as  exp(-Bt).  From  this  it  follows  that  the  relative  dielectric 
constant  is  then  continuous  everywhere.  Furthermore,  it  will  be  assumed 
that  the  magnetic  field  has  only  one  component  and  that  the  latter  is 
transverse  to  the  direction  of  propagation  which  has  been  chosen  to  be  x. 
Finally,  the  problem  is  transformed  into  a  two  dimensional  one  by  requir¬ 
ing  that  all  variations  along  the  z-axis  to  be  zero.  The  reader  is 
referred  to  Fig.  1,  p.22. 

Since  the  relevant  equations  for  the  regions  x<0  and  x>t  are 
special  cases  of  the  equations  describing  the  compressible  inhomogeneous 
plasma,  the  latter  region  will  be  considered  first. 

3.  DERIVATION  OF  THE  EQUATIONS  DESCRIBING  AN 
INHOMOGENEOUS  COMPRESSIBLE  PLASMA 

The  linearized  equations  describing  a  warm,  non-uniform,  source¬ 
less  plasma  as  formulated  by  Cohen  [1962)  are 
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where  Nq  and  v^  are  the  average  density  and  rms  velocity  of  the  electrons; 

— * 

m  and  e  are  the  mass  and  charge  of  the  electrons;  E  and  H  are  the  electro- 
magnetic  field  vectors;  and  v  and  are  the  velocity  vector  and  density 
of  the  electrons  imparted  by  the  fields.  In  addition  a  time  dependence 
of  the  form  e  iwt  has  been  assumed  and  suppressed  throughout. 

The  velocity  vector  appearing  in  (2)  may  be  eliminated  i>y  use  of 
(6)  yielding 
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Multiplying  (8)  by  Vx,  using  (1)  and  (3)  and  expanding  produces 
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Similarly  the  electric  field  vector  appearing  in  (2)  may  be 
eliminated  by  the  use  of  (6).  Multiplication  of  the  resultant  equation 
by  V*  produces 
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An  expression  for  V*(Nov)  may  be  obtained  by  solving  (7)  for  and 
substituting  the  resultant  expression  into  (5)  producing 
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The  last  term  appearing  in  (12)  is  a  static  term  and  is  henceforth 
dropped.  Upon  substituting  of  the  time  varying  portion  of  (12)  into 
(11)  and  expanding  one  obtains 
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An  expression  for  the  electric  field  may  be  obtained  by  solving 
(8)  for  E.  Then 
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Similarly,  an  expression  for  v  may  be  derived  by  elminating  £  from  (2) 
and  (6) .  Then 
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Hence  the  plasma  is  completely  specified  once  H  and  P  defined  by  (10) 
and  (13)  are  known. 


4.  SOLUTIONS  OF  THE  COUPLED  SYSTEM 
Let  us  expand  (10)  and  (13)  in  the  cartesian  coordinate  system 
for  the  case  where 
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In  general  the  y-variation  can  be  removed  by  a  Fourier  trans¬ 
form.  However,  since  we  are  dealing  with  the  propagation  of  plane  waves 
through  the  plasma  we  can  simply  assume  a  y-variation  of  the  form 
exp(iay)  where  a  «  kQ  sin  0q  and  0q  is  the  angle  of  incidence  measured 
from  the  normal  (Hessel  et  al  1962].  Equations  (17)  and  (lb)  become 
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where  ana  correspond  to  the  differential  operators  in  brackets. 
Equations  (19)  and  (20)  are  more  readily  solved  if  they  can  be  decoupled. 
Multiplying  (19)  by  produces 
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It  is  evident  tnat  P  may  be  eliminated  from  (21)  by  means  of  (20)  if  L, 
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Before  proceeding  with  the  solution  of  (23)  and  (24)  let  us  return  once 
more  to  (19)  and  (20)  to  see  if  an  additional  simplification  can  be 
effected.  For  the  exponential  profile  (22),  (19)  and  (20)  become 
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where  p  «  —  P  • 

An  equivalent  system  may  be  generated  by  the  addition  anc  subtraction  of 
(23)  and  (26)  which  has  tne  form 
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Unlike  Eqs .  (23)  and  (26),  (27)  and  (2b)  permit  a  comparison  between  the 
relative  magnitudes  of  kQ2  and  k^2  .  Let  us  consider  the  ratio 


P 


where  v  is  the  r.m.s.  velocity  of  electrons  and  c  the  velocity  of  light. 
P 

Normally  this  ratio  will  be  considerably  smaller  t.ian  unity,  then  to 
first  order,  any  contribution  due  to  kQ2  may  be  neglected.  The  effect  of 
neglecting  kQ2  in  (27)  and  (28)  is  equivalent  to  dropping  the  term 
involving  k  2  in  (23).  This  result  may  easily  be  verified  by  the  addi¬ 
tion  and  subtraction  of  (27)  and  v28)  and  the  subsequent  use  of  (29)  tnus 
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recovering  our  original  system  (25)  and  (26)  with  the  term  involving 
kQ  absent.  Thus  we  are  faced  with  the  task  of  solving 


and 


2 

A  l  +  B  dx  ■  * 
dx 


H  -  —  aUP 
2 


(31) 


^2  *  2  -Bx  2 


idx 


~  +  B  ” —  +  k  e 
2  dx  p 


P  »  ^  aBH 
e  z 


(32) 


By  means  of  the  transformation 


P  -  e“B^  F(x) 


03) 


„  -Bx  , /  % 
H  -  e  —  (.(x) 


(3M 


(31)  and  (32)  become 


$  ~  <’2  ♦  t)0 


wtn 


oBF 


(35) 


and 


d“F  .  ,,  2  -Bx  2  B2. 

— 7  +  <k  e  "  a  ~  T  )F 
dx2  P  4 


aim 


aBG 


(35) 


Lquations  (35)  and  (36)  may  be  decoupled  according  to  the  method  indicated 
oy  (23)  and  (>4).  Then 
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Before  applying  the  method  of  Frobenius  to  each  of  the-  otcouplcd 
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equations  it  is  convenient  to  remove  the  factor  e  .  To  this  enu  let 
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Assuming  power  series  solutions  of  the  form 
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one  obtains  upon  substitution  into  (40)  and  (41)  the  following  recursion 
relationsnips 
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where  s  is  determined  from  tne  indicial  equation 


(43) 
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Since  two  of  the  roots  of  the  indicial  equation  are  also  zeros  of  the 
denominators  of  the  recursion  relations  for  the  case  n  »  2,  both 
solutions  of  the  first  and  second  kina  must  be  constructed.  As  this 
procedure  is  rather  well-known,  only  the  final  results  will  be  pre¬ 
sented.  Thus  the  solutions  generated  by  the  method  of  Forbenius  are 
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The  above  solutions  hold  for  the  case  where  iLSl  is  not  equal  to  an 
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integer  or  zero. 

Substituting  the  generated  series  into  (31)  and  (32)  reveals 
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B  aim  B 


com 


however,  may  be  found  by  forming  linear  superpositions  of  the  various 
series  generated.  Thus,  if 


H_(x)  =  ~  y  H.(x)  +  y_  H_(x) 


am 


'A  A 


B  c 


P,  (x)  «  1  P  (x)  +  Y  P  (x) 
3  c  A  U  c 


(46) 

(47) 


then  a  substitution  of  (46)  and  (47)  into  (31)  and  (32)  reveals  that 


*  y 


(B  -  2a) 


c  A  2aB(B  +  2:0  8 

YD  =  '  yb. 

Hence,  the  various  y’s  cannot  be  determined  uniquely. 
(48)  and  (49)  into  (46)  and  (47)  produces 


(48) 


(49) 


Substitution  of 


P3  ~  YB 


kp  <B-2a)  ^ 
2aB(B+2a)  ?A  "  ?C 


+  yapa 


—  Yj,H  +  -£-  y ,H 

Be  A  A 


um 


(50) 

(51) 


or  alternately 


P3  “  “  YB  Pc  +  Yc  PA 


k  (B-2a) 

H,  «  —  Y„  H  -  ^TTTTT-T  H  +~yH, 
3  atm  Bl  c  2aB(B+2a)  AJ  com  c  A 


The  last  terms  in  (50)  to  (53)  may  be  dropped,  since  neing  already 
solutions  of  the  system, they  add  nothing  new.  That  tne  alternate  forms 
are  equivalent  and  that  eitner  representation  may  be  used  for  too  prob¬ 
lem  at  hand  will  be  demonstrated  shortly.  Similarly,  if 


n,(x)  “  “  *fr  Y,  H  (x)  -  ~  y,  b  (x) 
4  u*n  G  B  ojm  1;  b 


p4(*>  -  rt  PB  +  yF  P„ 


k  (B  +2a) 

y  n  Y  +  *  . . —  Y 

’t  tG  -2aB(B-2u) 


YH  "  "  YF 


k  (B  +  2a) 

P  w  Y  —2 -  P  -  P  +  Y  P 

4  YH(-2aB(B-2a)  B  ol  G  B 


Y  H  -  —  y  H 
rH  com  D  com  'B  B 


Y  P  +  Y  P 
TH  D  't  8 


,'AUl 
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..  k  2CB  +  2a) 

H  m  ~  P  -  -  - .  -  -  u  _  ll 

4  cum  -2aB(B-2ct)  B  U 

. 

Again  the  last  terms  appearing  in  (59)  through  (61  may  be  dropped  since 
they  are  already  solutions  of  the  coupled  system. 

The  sets  of  solutions  satisfying  (31)  and  (32)  are  summarized 

below 


-  y_  ■— 

E  wm  B 


(61 


Pl<*> 

H1(X) 

P2(x) 

U2(x) 

P3(x) 

H3(x) 


“  PA(x) 


m 


e 

cum 


V*> 


■  V*> 


-  -  rr  H  <x> 
cum  o 

k  2(B-2c) 

‘  2&><Th*.7  PA(X)  -  Pc(x) 


e 


H  (x) 
c 


k  2(B+2a) 

Vx)  "  -2aB(B-2a)Vx)  "  PD(x) 

VX>  "  "  £  HD(X)  . 


(62) 


The  complete  solution  describing  the  pressure  and  magnetic  field  for  the 
region  0  £  x  £  t  is  obtained  by  forming  a  linear  superposition  of  all  the 
solutions  tabulated  above.  Thus 


H(x)  -  C1H1(x)  +  C2H2(x)  +  C3H3(x)  +  C4H4(x) 

and 

P(x)  -  C1P1(x)  +  C1P1(x)  +  C^P^(x)  +  ^(x) 


(63) 
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t 


Recalling  that  the  C’s  are  yet  unspecified  the  substitution 

(65) 

(66) 

transforms  (65)  and  (66)  into  the  form  which  would  have  resulted  if  the 

alternate  forms  for  P~,  P.  and  H.  had  been  used.  Thus  it  can  be 

3  4,  3  4 

concluded  that  the  various  forms  are  equivalent. 

5.  THE  BOUNDARY  CONDITIONS 

Having  obtained  solutions  to  the  coupled  set  of  electromagnetic 
and  hydrodynamic  equations  a  consideration  of  the  boundary  conditions  at 
the  interfaces  located  at  x  ■  0  and  x  ■  t  is  in  order.  For  the  electro¬ 
magnetic  fields  we  shall  require  that  the  tangential  components  of  the 
electric  and  magnetic  fields  are  continuous  at  the  interfaces.  In  order 
to  determine  the  behavior  of  P  at  the  free  space-plasma  interface  Eq.  (6) 
is  expanded  in  component  form.  Thus 

ap 

•r*N  e  E  -  iojtn  N  v  (67) 

3y  o  y  o  y 

and 

-  ■“  *  N  e  E  -  iwm  N  v  (68) 

dx  o  x  o  x 

2 

where  N  *  tonUji  (1  -  e  Bx)  , 
o  r 

e 

Since  the  electron  density  is  equal  to  zero  for  x  equal  to  zero  ano 
the  components  of  tne  electric  field  and  ion  velocity  are  expected  to  be 


k  (B-2o) 


2aB(B+2a) 


k  (B+2a) 

>•  «.  _£ — — _ - 

'2  -2aB(B-2a) 


I 

i 


I 
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bounded  it  follows  that 


2P 
3  x 

x»o+ 


0 


(69) 


r  “0  (70) 

x*o+ 

In  addition  we  shall  require  that  P  be  continuous  at  x  »  t. 


6,  THE  DETERMINATION  OF  THE  REFLECTION  AND  TRANSMISSION  COEFFICIENTS 
The  partial  differential  equation  describing  tue  magnetic  field 
in  free  space  is  obtained  from  (15)  for  the  case  e  *  1.  The  solutions 
are  well-known  and  the  expression  for  tne  magnetic  field  may  be  written 
as 

ik  cos©  x  -ik  cos©  x  ik  sin©  y 
rlzi  «  (e  +  Re  )  e  (71) 

where  ft  is  tne  reflection  coefficient  sought  and  ©^  the  angle  of  incidence 
measured  from  the  normal.  Using  (15)  it  follows  tnat  the  expression  for 
tne  tangential  component  of  the  electric  field  is 

k  cos©  ik  cos©  x  -ik  cos©  x  ik  sin©  y 

„  OO,  OO^OO^  OO'  /-,„x 

E  .  *  — - (e  -  Re  )  e  .  (72) 

yl  we 

o 


For  the  region  occupied  by  the  compressible  inhomogeneous  plasma  we  have 


Hz2  "  C!H1(X)  +  C2H2<X)  +  C3H3(X)  +  Wx) 


ik  sin©  y 

0  °  /71\ 

e  (73) 


P  - 


ClPl(x)  +  C2P2(x)  +  C3?3(x)  +  C4P4(x)'  6 


ik  sin©  y 
o  o 


(74) 


The  tangential  component  of  the  electric  field  E  £  ke  calculated 


1? 


from  (15).  By  making  use  of  (73)  and  (74)  this  may  be  expressed  as 

-1 


~y2  weoe(x) 


c.  (h!  -  —  p.)  +  c,(h;  -  —  p9)  +  c,(n;  -  —  P  , 

11  ura  1  22  wm  2  33  um  3 


+  c, (h:  -  p, ) 

4  4  10m  4 


(73) 


where  the  prime  denotes  differentiation  with  respect  to  x  anc  a  *  kQsin©o. 

The  equations  describing  the  pressure  and  the  magnetic  field  for 
tne  region  occupied  by  the  Homogeneous  compressible  plasma  may  be  ob¬ 
tained  from  (17)  and  (18)  setting  e  *  ■  exp(-Bt)  »  constant. 


Then 


d?'H 


dx 


f  +  (ko  E3  '  “Xs  "  ° 


(70 


2 

a  P 


dx 


T  *  <kp2  e3  -  ‘X  ■  °  • 


(77) 


Since  the  region  for  x  >  t  is  semi-infinite  only  outgoing  waves  will  be 
present.  Tnus  solutions  for  tnis  region  arc  of  tne  form 


u  _  +i -  a  (x-t)  ik  sin©  y 
H  *  Te  o  3  eo  o'' 
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(78) 


n  _  2c.  -  a2  (x-t)  ik  sin©  y 

«  Qe  p  3  e  o  o/ 


CO 


As  before  the  tangential  component  of  the  electric  field  obtainec 

from  (15).  Making  use  of  (78)  and  (79) 


>3 


~1 


oe3  a  T  ii/C2e„  -a2  (x-t) 
Te  o  3 


{  ueoe3 


+  -J£SL  , 

2 

io  me  e. 
o  3 


i/k  2e  -  a2  (x-t) 
Qe  p  3 


ik  sin©  y 
e  o  o' 


(80) 
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Application  of  the  indicated  boundary  conditions  yields  (81)  where  the 
desired  coefficients  may  be  determined  by  applying  Cramers  rule. 

7 .  NUMERICAL  RESULTS 

For  purposes  of  numerical  evaluation  the  reflection  and  trans¬ 
mission  coefficients  may  be  expressed  as 

W,  -  k  cos©  >4  -  a*’  W.  -  i(k  cos©  W_  +  »4^e_  -  k'  ) 

„  looo3  4  o  o  3  o  3  2 


i0X~2'  2 


W,  +  k  cos©  /k  -  a*-  W,  +  i(k  cos©  W_  -  /k^e0  -  W  ) 
1  o  o  o  3  4  c  o  3  o  3  2 


(62) 


i2k  cos©  We 
o  o  5 


W,  +  k  cos  <4  "2c.  -  a"  W.  +  i(k  cos© 

1  o  oo3  4  o  o3 


W„  -  A  2e,  -  a2  WJ 


o  3 


2 

(83) 


where 


H’(o) 

H’(o) 

«2  (o) 

H;(o) 

0 

Px(o) 

P2(o) 

P  3(o) 

P4(o) 

0 

p;(o) 

P'(o) 

P^(o) 

PA(o) 

U 

Px(t) 

P2(t) 

P3(o) 

P4(t) 

-1 

Hx(t) 

-  — p,(t) 
aim  1 

H'(t) 

"  2k<*>  «3(t) 

ea„ 

- P. 

uxn  . 

s(t)  H^(t) 

-  £2p (t)  — 
cum  4  mm 

(84) 

H[(o) 

H'(o) 

H’(o) 

H^(o) 

w  - 

Px(o) 

P2(o) 

P3(o) 

P4(o) 

(85) 

z 

P[(o) 

P^(O) 

P^(o) 

?l(o) 

Hx(t) 

H2(t) 

H3(t) 

J 


8,(0) 

& 

N) 

O 

H3(o) 

v°> 

0 

PL(o) 

P2(o) 

P3(o) 

P4C0) 

0 

w  » 

3 

P^(o) 

p;(o> 

P3(o) 

*;<•> 

0 

P1(t) 

p2(t) 

P3(t) 

-1 

H2(t)“  fj“Mt) 

i.  wm  / 

H' (t)-  ~P^(t) 
3  M  3 

KM-  SV** 

eo 

wm 

(86) 


1^(0) 

«2  (°) 

h3(°) 

H4(o) 

Pl(o) 

P2(o) 

P3(o) 

P4(o) 

V 

P|(o) 

P](o) 

P}(o) 

Pi(0) 

H^t) 

IL(t) 

H3(t) 

H4(t) 

P^o) 

P2(o) 

P3(o) 

P,(o) 

0 

P|(o) 

P2(o) 

P3(o) 

PJ<  0) 

U 

»L(t) 

H2(t) 

h3(t) 

h4(0 

0 

PL(t) 

P2(t) 

P3(t) 

P4(t) 

-1 

&!<*> 

H~ (t)-  ”P7(t) 

2  wm  2 

H*  (t)—  ”P,(t) 
3  um  3 

H!  (t)-  (t) 

4  wm  4 

ea 

wm 

(88) 

It  should  be  noted  that  all  of  the  entries  In  tne  above  determinants  are 

real.  Furtnermore,  expressions  for  tne  reflection  and  transmission  coef- 

2  2 

ficients  for  tne  case  where  3  -  a  <  0  may  be  obtained  by  replacing 

the  factor  (k^  appearing  in  (62)  and  (83)  by  +i(a^  -  k^tj)^^ 

While  all  of  tne  series  solutions  converge,  their  rate  of  conver- 

2 

gence  depends  critically  upon  the  magnitude  of  the  factor  (k^/B)  .  Using 
(14)  and  the  expression  for  the  acoustic  velocity  in  an  electron  gas  at  a 
temperature  T 


3kT 


(89) 
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1 


t 


it  follows  that 


4ir2m  f2 
3k  TB2 


(90) 


where  m  is  the  electronic  mass,  k  the  Boltzmann  constant,  and  f  the 
frequency  of  the  incident  electromagnetic  radiation.  For  a  frequency  of 

2 

1GHz,  a  temperature  of  2000*K,  and  B  equal  to  unity,  tne  value  of  (k^/B) 
will  be  of  the  order  108.  Consequently,  all  of  the  series  will  converge 
rather  slowly  ana  require  the  calculation  of  many  terms.  Furthermore, 
these  terms  rapidly  reach  magnitudes  which  exceea  the  capabilities  of 
most  computers.  While  it  is  believed  that  this  problem  can  be  overcome 
by  special  programs,  a  lack  of  time  prevented  any  efforts  in  this  direction. 

In  order  to  obtain  some  idea  of  the  behavior  of  the  functions  in- 
2 

volved,  (kp/B)  was  set  equal  to  unity.  Then  for  a  frequency  of  1GHz, 

B*50,  and  the  thickness  of  the  compressible  inhomogeneous  plasma  layer 
equal  to  .014  meters,  a  reflection  coefficient  of  about  unity  and  a  trans¬ 
mission  coefficient  effectively  zero  were  observed  for  all  angles  of  inci¬ 
dence  between  5  and  85  degrees.  It  must  be  pointed  out  that  the  assumed 
numerical  values  actually  describe  a  plasma  at  a  temperature  so  high  that 
it  is  no  longer  consistent  with  the  governing  equations  and  the  assumptions 


behind  them. 
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